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Abstract 

We address the thermal history of the Earth after the Moon-forming impact, 
taking tidal heating and thermal blanketing by the atmosphere into account. The 
atmosphere sets an upper bound of ~ 100 W/m^ on how quickly the Earth can cool. 
The liquid magma ocean cools over 2-10 Myrs, with longer times corresponding to 
high angular-momentum events. Tidal heating is focused mostly in mantle mate¬ 
rials that are just beginning to freeze. The atmosphere’s control over cooling sets 
up a negative feedback between viscosity-dependent tidal heating and temperature- 
dependent viscosity of the magma ocean. While the feedback holds, evolution of the 
Moon’s orbit is limited by the modest radiative cooling rate of Earth’s atmosphere. 
Orbital evolution is orders of magnitude slower than in conventional constant Q mod¬ 
els, which promotes capture by resonances. The evection resonance is encountered 
early, when the Earth is molten. Capture by the evection resonance appears certain 
but unlikely to generate much eccentricity because it is encountered early when the 
Earth is molten and Qe 3> Qd • Tidal dissipation in the Earth becomes more effi¬ 
cient (Qe Q(i) later when the Moon is between ~ 2Oi?0 and ~4Oi?0. If lunar 
eccentricity grew great, this was when it did so, perhaps setting the table for some 
other process to leave its mark on the inclination of the Moon. 
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Highlights 

• We address the thermal history of the Earth after the Moon-forming impact. 

• Thermal blanketing by Earth’s atmosphere limits cooling to ~ 100 W/m^. 

• Feedback between tidal heating and viscosity tethers the Moon to Earth. 

• Slow orbital evolution promotes capture by resonances. 

• Tidal Q’s evolve differently in Earth and Moon 
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1. Introduction 


The leading explanation for the origin of the Earth-Moon system is that it formed 


when two planets collided 

(Hartmann et al. 

1986 

Benz et al. 


Stevenson 

1987 

Canup and Righter 

2000 

Canup 

2004 

Reufer et al., j2012 ; 

Cuk and Stewart 

2012 


Canup, 2012). The event is widely known as the Moon-forming impact. But it might 


have been called the Earth-forming impact. Indeed, in some scenarios (e.g., Canup 


2012), neither of the two colliding planets is markedly more a proto-Earth than the 


other, in which case the Moon-forming impact is literally the Earth-forming impact 
as well. 

Here we accept as an initial condition that a Moon-forming impact took place. 


We focus on the angular-momentum-conserving canonical model (Canup, 2004), al¬ 


though we also consider some aspects of recent non-angular-momentum-conserving 


alternatives (Reufer et al 


2012 


Cuk and Stewart 


2012 


Canup, 2012) that may 


better address isotopic evidence that the Earth and Moon share the same source 


e.g., Pahlevan et al., 2007 Zhang et ah, 2012). 


Our subject here is the earliest thermal evolution of the Earth after the impact. 
Earth’s cooling was greatly slowed by the thermal blanketing effect of its atmosphere, 
which was thick and deep and very opaque, and Earth’s interior was strongly heated 
by tidal flexing, because the Moon was close by and the Earth was spinning rapidly. 
We stress here that tidal heating and thermal blanketing were linked, because radia¬ 
tive cooling by the atmosphere limits how big tidal heating can be. Hence we are 
also concerned here with the early tidal evolution of the Moon’s orbit, which we will 
treat as simply as possible. We hrst describe the thermal model of the cooling Earth 
in Section We then describe tidal heating and tidal evolution in Section 


2. The cooling Earth 

2.1. Initial conditions 

A reasonable initial condition on the Earth after the Moon-forming impact is that 


it behaved like a hot liquid (Zahnle et al., 2006; Elkins-Tanton, 2008). The impact 
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itself melted and partially vaporized the side of the Earth that was directly struck. 


But even if half the mantle survived the shock of collision unmelted (Tonks and 


Melosh, 1993), it would have collapsed under its own weight on a Rayleigh-Taylor 


instability timescale of the order of 


r 


pp 


( 1 ) 


gApL' 

where L is the relevant length scale, g gravity, Ap/p the density contrast between 
the solid and the melt, and u the kinematic viscosity of the solid. For reasonable 
choices {Ap/p = 0.02, L = 3 x 10® cm, u = 10^® cm^/s, the latter appropriate to a 
warm silicate planet during accretion), the solid hemisphere would collapse in just 
5 years and in the process releasing more than enough gravitational energy to melt 
what had not yet melted (Tonks and Melosh, 1993). Other factors — the impactor’s 
shock-heated molten iron plunging through the mantle to join with Earth’s core and 
the superheated state of the upper core immediately thereafter — also suggest that 
a liquid planet is a plausible initial condition. 

It is important here to be careful dehning what is meant by melting. Melting of a 
silicate is characterized by three temperatures: a liquidus Tuq hotter than which the 
silicate is fully molten; a solidus Tgoi colder than which the silicate is fully solid; and 
a critical temperature Tent in between at which the rheology of the material changes 
from that of a solid with melt percolating through it (T < Tent) to a liquid carrying 


suspended solids (Abe, 1997 Solomatov, 2009). The critical temperature Tent marks 
a phase transition in viscosity, which best corresponds to what is colloquially meant 
by a melting point. When the mantle collapses under its own weight, it is viscosity 
that converts that motion into heat. Viscous dissipation can only raise the temper¬ 
ature of the mantle to Tent- To raise the temperature of such a lower mantle to the 
liquidus requires another mechanism, perhaps a prompt effect like the impactor’s 
shock-heated molten iron dispersing in the mantle before it joins with Earth’s core, 
or perhaps a delayed effect stemming from the mantle being strongly heated from 
below by the core. 
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Although very important to geochemistry, for our purposes it is not very impor¬ 
tant whether the mantle be fully or only partially molten, provided only that initially 
it behave like a low viscosity liquid. We therefore begin our study with a liquid man¬ 
tle and a surface temperature of 3000-4000 K some 100-1000 years after the impact, 
by which point we can hope that early transients have settled down. On the other 
hand, these transients may last long enough to influence moon formation, which is 


also generally regarded as very fast (Canup, 2004 Salmon and Canup, 2012). 

A second initial condition is the presence of a substantial atmosphere. Earth 
today has about 270 bars of H 2 O and 50 bars of CO 2 near the surface, the former 
mostly in oceans and the latter mostly in carbonate rocks. Similarly large reservoirs 


are held in the mantle (Sleep et al., 2001). Thus it is reasonable to start with 100- 
1000 bars of H 2 O and CO 2 in the atmosphere. The water vapor and carbon dioxide 
would be supplemented by smaller amounts of CO, H 2 , N 2 , various sulfur-containing 


gases, and by a suite of geochemical volatiles evaporated from the magma (Schaefer 


et al., 2012; Fegley and Schaefer, 2013). 


A third initial condition is the existence of the Moon. We start the Moon with its 
current mass and density at the relevant Roche limit of 2.9i?0. Although accretion 
theory suggests that the Moon could accumulate in less than a year, the need to ra¬ 
diate the energy liberated by lunar accretion likely prolonged the process by decades 


or even centuries(Thompson and Stevenson, 1988 Ida et al., 1997 Genda and Abe 


2000; Canup, 2004 Salmon and Canup, 2012). During this time interactions between 


the accreting Moon and the unaccreted debris disk may raise the starting point of 


the Moon beyond the Roche limit to perhaps ~ (Salmon and Canup, 2012). For 
our purposes starting the Moon nearer the Earth is more interesting because there 
is more tidal heating and more resonances to encounter as the Moon’s orbit evolves. 

A fourth initial condition is the initial angular momentum of the Earth-Moon 
system. In canonical models this is held constant, so that with the Moon forming 
at 2.9i?0, the Earth starts with a 5.0-hour day. Some recent Moon-forming models 


posit loss of angular momentum from the Earth-Moon system (Cuk and Stewart 
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2012; Canup, 2012). These begin the Earth with a ~ 3-hour day, and thus with 


signihcantly more rotational energy available to be dissipated as tidal heating. 


2.2. Interior Cooling 

The total thermal energy E of the Earth is the sum of the thermal energy in the 
interior Ei and the thermal energy in the atmosphere Eatm- Our interest is almost 
entirely in the time rate of change of E, 

E = —-|- ttR^^Eq (1 — A) -1- Eude- (2) 

Equation describes the Earth’s cooling rate as the difference between thermal in¬ 
frared radiation emitted to space by the atmosphere and the sum of net insolation 
(radiative heating), tidal heating, and heating from the decay of radioactive elements. 
Thermal radiation (radiative cooling) to space is described by an effective radiating 
temperature which is determined by the thermal blanketing effect of the at¬ 
mosphere and the surface temperature T^; a is the Stefan-Boltzmann constant and 
is the radius of the Earth. The relation between and Tg — the greenhouse 
effect — is discussed in the next section. Radiative heating is described by the solar 
constant Fq ~ 1000 W m ^ ca. 4.4 Ga and an albedo A. We leave A = 0.3 as today. 
Tidal heating Etide is a big term addressed in Section below. With radioactive 
heating ignored, the geothermal heat flow E is equivalent to the rate of cooling of 
the interior, 

AttRIE = Etide - Ei. (3) 


2.3. Radiative Cooling 

It is well known that an atmosphere’s thermal blanketing effect prevents a magma 
oceauQfrom cooling rapidly. A string of increasingly sophisticated models, beginning 
with|Matsui and Abe (1986) and proceeding through Abe and Matsui (1988); Zahnle 


^The term “magma ocean” is broadly used in the literature to describe a mantle with significant 
partial melt. 

























et al. 

(1988) 

, Abe 

(1997) 

Zahnle et al. 

(2006) 

Elkins-Tanton 

I2(l(ls 

2012) 

Hamano 

et al. 

(2013) 

, and 

Lebrnn et al. 

(2013) 

have addressed some of the conseqnences of 


thermal blanketing by steam atmospheres over the Earth’s cooling magma oceans. 
Mnch of the progress made in these stndies has been in constraining how the magma 
ocean freezes, how long it takes to freeze, and how, when, and what it degasses. All 


of these models presnme H 2 O-CO 2 atmospheres, and all bnt Zahnle et ah (2006) 
ignore the Moon. 


A different set of stndies has addressed the composition of the atmosphere (Schae¬ 


fer and Fegley 2010 Schaefer et ah 2012 Fegley and Schaefer , 20131 Lnpn et ah 


2014). Shortly after the impact the atmosphere over the lava is very hot and con¬ 


tains a great many geochemical volatiles that can evaporate from a magma ocean, 
snch as snlfnr, sodinm, and chlorine, and many other componnds of H, C, N, and 


O in addition to H 2 O and CO 2 (Schaefer et ah, 2012 Fegley and Schaefer, 2013). 
The geochemically enriched atmosphere is more opaqne than jnst H 2 O and CO 2 at 
magma temperatnres and thns the planet cools more slowly and spends more time 
with a lava snrface than previons models wonld predict. 

To address these issnes we make nse of new 1-D non-gray radiative-convective 
atmospheric strnctnre models developed to compnte the Inminosity and spectral 


characteristics of Earth-like planets after giant impacts (Lnpn et ah, 2014). As the 


Lnpn et al. (2014) models are discnssed fnlly elsewhere we will not provide fnll doc- 


nmentation here, bnt we do note some key featnres. The new models inclnde many 
of the molecnlar and atomic opacity sonrces that wonld be abnndant in eqnilibrinm 
with a hot rock or magma snrface with the composition of the bnlk silicate Earth 


(BSE, Lnpn et ah, 2014). The thermal strnctnre models are compnted as eqnilib- 
rinm models in which thermal radiation emitted by the atmosphere is balanced by 
the snm of snnlight absorbed and geothermal heat extracted from the cooling Earth. 
This assnmption is very good becanse the heat capacity of the atmosphere is very 
small compared to the heat capacity of the planet. A more qnestionable assnmption 
of the model is that it keeps the atmosphere in chemical eqnilibrinm with the bnlk 


7 










































































































silicate Earth at any temperature — i.e., it assumes that CH4 and NH3 form when 
equilibrium chemistry would favor them. The assumption is reasonably good at high 
temperatures (at which CH4 and NH3 are rare in any case), but for temperatures 
much below 1000 K equilibrium would not be expected unless effective catalysts 
were present among the cloud particles. At low temperatures, equilibrium predicts 
that CH4 and NH3 would be present and effective as greenhouse gases, and thus 


the published Lupu et ah (2014) model may overestimate the greenhouse effect and 
underestimate Teg when surface temperatures are below 500 K (above which H 2 O is 
abundant). On the other hand, we are missing many opacity sources that would be 
present at higher temperatures, and thus in all likelihood the model underestimates 
the greenhouse effect when surface temperatures are above 1000 K. In practice we 
address the consequences of uncertain opacities by varying the total pressure of the 
atmosphere. 

New calculations of Teg for cloud-free runaway greenhouse H 2 O-CO 2 -N 2 atmo¬ 


spheres (Goldblatt et ah, 2013; Kopperapu et ah, 2013) hud that Tes = 265 during 
the runaway greenhouse phase of planetary cooling We expect that a deep CO 2 atmo¬ 
sphere will remain after the ocean of water condenses. This state was not addressed 
explicitly by Lupu et al.| (2014). We therefore retrohtted a smooth transition to the 
~500 K surface temperature expected under 100 bars of CO 2 (Zahnle et ah, 2006). 
What happens to the CO 2 at later times on longer timescales is discussed elsewhere 


(Sleep et ah, 2001, 2014). 


Physical parameters 

For present purposes a ruthlessly simplihed model of the Earth is the best choice. 
We characterize the Earth by a surface temperature Tg and an interior temperature 
Tj. In effect T* is the potential temperature, the temperature that every parcel in an 
adiabatic mantle would have were it brought to the surface. We also assume that 
the mantle is uniform, in the sense that it has the same viscosity at all depths and 
that it freezes everywhere at the same time. This assumption is roughly equivalent 
to assuming that the adiabat and the melting curve are parallel. We take Tuq = 1800 






















K and Tso\ = 1400 K. These are arbitrary values; our results are insensitive to them. 
We assume that the melt fraction (j) is linear in Tj between Tuq and Tsoi, 

0=1 Ti> Tliq 

0 = {Ti — Tsoi) / (Tliq — Tsoi) Tliq > Ti > Tgol 

0=0 Tsoi > Ti. (4) 

For Ti > Tliq and Ti < Tsoi, the heat capacity of the Earth is 

= C'^,manll^man + Cy,cor^cor = 6.2 X 10^*^ J/kg/K, (5) 

in which the masses of the Earth, the mantle, and the core are M^, and Moor, 

respectively. The heat capacity of silicate is approximated by Ct,,man = 1200 J/kg/K 
and that of iron by Cv,cor = 650 J/kg/K. The iron core does not freeze on the time 
scale of interest. While silicates freeze, the effective heat capacity of the Earth is 

= C„®A4 + = 1.0 X 10“ J/kg/K (6) 

where the heat of fusion is Qm = 4 x 10® J/kg. With the above simplifications, the 
relation between Ei and Ti is 

Ei = Cy^M^Ti 
Ei = C^^M^Ti 

Including the atmosphere, the total E is 

E = Ei + Ctj^atmlUatmTs, (8) 

in which Tg is related to T* by the equations of parameterized convection described 
in the next section. The thermal energy in the atmosphere is a small term. 


Ti > Tliq, ^sol > Ti 

Tliq > Ti > Tsoi. (7) 


2.5. Parameterized convection 

We follow previous workers by using parameterized convection to link the cooling 


rate to the temperature and heat generation inside the Earth (Solomatov, 2009 
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Lebrun et al., 2013). In parameterized convection, the geothermal heat flow F 

kc {Ti 


F = C. 


T, 




(9) 


is related to the size and viscosity of the interior and the temperature gradient across 
the surface boundary layer through the Rayleigh number Ra, 

a^g {Ti - Ts) 


Ra = 


KV 


( 10 ) 


Equations and 10 can be combined to give a closed-form expression for F{Ts,Ti), 


F = Cok, (T,-T,)'+”L3n-i 


O.'uQ 

Kiy{Ti) 


( 11 ) 


We use kc = 500 J/m^/s/K for the thermal conductivity; is the volume thermal 
expansivity, g the gravity, and n is the thermal diffusivity (k = kc/C^^caan/p)- We 
will use = 2.4 x 10“^ ^ = 10“® m^/s, and 5 ^ = 10 m/s^. The characteristic 

distance scale L is usually identified with the thickness of the mantle. In the soft 


turbulent regime we take Cq = 0.089 and n = 1/3 (Solomatov, 2009). With n = 1/3, 
L cancels out, so that the heat flow is determined only by the properties of the 
boundary layer. 

The kinematic viscosity u (m^/s) is a very strong function of temperature R 
and therefore of central importance to convection and tidal heating. We follow 


Abe (1997); Solomatov (2009); Lebrun et al. (2013) by treating the temperature 


dependence of viscosity in two regimes, the hotter resembling a liquid and the cooler 
resembling a solid, the regimes divided by Tent. The phase transition is a sharp 


function of the fraction 0 of silicate that is liquid. In the liquid, Solomatov (2009) 
parameterizes viscosity by 

2.5 


u {Ti > Tcrit) = 6 X 10 ® exp 


4600 K 


1 - 


-'crit 


mVs. ( 12 ) 


Ti — 1000 Ky \(y — (ycrit 
This expression implies that z/ —)■ cx) as 0 —>■ 0crit; he., it is singular at 0 = 0crit- In 
the solid, 

Tsoi — Ti 


z/ (Tj <Tcrit) = 1 X 10 ^'" exp 


58 K 


exp 


-a 


m^/s. 


(13) 
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with = 26 (Solomatov, 2009). Abe (1997) places the boundary at 0 = 0crit = 
0.4. For 0crit = 0.4, Tent = 1560 K. The particular values of 0crit and Tcrit are 
unimportant. Figure [T] illustrates the situation. 



Temperature 


Figure 1: Semi-quantitative cartoon illustrating the rheological transition in a cooling magma 
ocean. Cooling proceeds from left to right. At the liquidus refractory crystals begin to freeze 
out. With further cooling the crystal fraction increases and the liquid fraction 0 decreases. At the 
rheological transition 0 = 0c ~ 0.4. At the solidus the melt fraction 0 = 0. The corresponding 
kinematic viscosity v is shown on a logarithmic scale. The discontinuity in z/ at 0 = 0c (T = Trit) 
corresponds to a melting point. 


2.6. Moonless Earths 

Impacts on the scale of the Moon-forming impact need not create a moon. It is 
illustrative to consider a case without the Moon for comparison with previous work 
(Hamano et ah, 2013; [Lebrun et ah , 2013). Computed temperatures and heat flow 
histories for a 100 bar BSE-equilibrated atmosphere are shown in Figure In these 
simulations the initial temperature is set at Tj = 3500 K at t = 100 years. 

Cooling takes place in two sharply dehned regimes (Figure]^. In the hrst regime 
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Figure 2: Thermal evolution of an Earth after a non-moon-forming giant impact of the same 
energy as the canonical Moon-forming impact. Results are shown for a 100 bar BSE atmosphere. 
Internal temperature Ti (mantle), surface temperature Tg, and effective radiating temperature Teff 
are plotted against the left-hand axis. The liquidus and solidus are indicated. Geothermal heat 
flow F (dotted curve) is plotted against the right-hand axis. The heat flow plateaus at F ~ 110 
W/m^ while the Earth passes through the runaway greenhouse, a phase that lasts from ^0.5 Myr 
to ~ 1.7 Myr. This period corresponds to a liquid mantle. It ends when the mantle becomes too 
viscous to sustain F ~ 110 W/m^, at which point Tg drops to low values, here ~500 K under ~ 100 
bars of CO 2 . The mantle remains partially molten for ^ 100 Myrs. The latter is what is often 
meant by “magma ocean.” 


the three leading terms in the planet’s energy budget, 

trT'j = (1 - y4) ^ 0/4 + F, (14) 

are of comparable magnitude, with the geothermal heat flow F ~ 110 W/m^ making 
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up the difference between net insolation ~ 170 W/m^ and radiative cooling at the 
runaway greenhouse threshold ~280 W/m^. In this regime the temperature contrast 
AT = Ti — Tg across the surface thermal boundary layer is small compared to Ti. 
This lasts about 1 million years. The transition to the second regime takes place 
when the mantle’s viscosity grows too great to sustain F ~ 110 W/m^, at which 
point the mantle undergoes the phase transition depicted in Figure After the 
phase transition heat flow is negligible (here ~ 20 W/m^ at first, but falling off fairly 
quickly) and the surface temperature is determined by insolation and the greenhouse 
effect of ~100 bars of CO 2 . 

As Matsui and Abe (1986) and Hamano et ah (2013) point out, the value F ^ 110 


W/m^ is particular to ancient Earth with insolation Fq = 1.0 x 10^ W/m^ and 
albedo A = 0.3. At Venus’s distance from the Sun, the heat flow required to support 
a runaway greenhouse atmosphere with A = 0.4 ±0.1 is of the order of T = 10 ± 40 
W/m^. Venus after the impact could have been held in a state with a molten surface 
for tens of millions of years or more, and if its albedo were less than A k. 0.38 it need 


not have emerged from this state until after its water was lost to space (Hamano et 


3. The Moon and Tides 


Tidal heating by the nearby Moon was a major term in early Earth’s energy 
budget. The heat comes from braking the rotation of the Earth. Tidal dissipation 
is usually described by the Q parameter, defined as the fraction of energy in the 


wave that is damped out during a single period (Goldreich and Soter, 1966). The 


value of Q can often be deduced empirically while presupposing nothing about the 
underlying physics of tidal dissipation. We will use Q, but for the congealing man¬ 
tle we will assume that tidal dissipation can be described by viscosity. This is a 
reasonable approximation as the mantle starts to freeze because ordinary viscous 


dissipation would be big (Moore, 2003). Where viscous dissipation dominates, Q 


can be computed (i.e., instead of assuming Q directly, we assume a particular physi- 
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cal mechanism for dissipation and a particular description of how viscosity evolves). 
But at very early times when the mantle is effectively inviscid, or at modern times 
when the mantle is solid, ordinary viscous dissipation is probably negligible and some 
other source of dissipation will likely be more important. In a hot molten Earth we 
use observed values of Q in fluid planets as a guide. For modern times we assume 
an asymptotic viscosity that brings the Earth-Moon system to its current state in 
4.5 Ga, it being understood that the bulk of tidal dissipation on the Earth today 
actually takes place at the interface between the solid surface and the liquid water 
oceans 


(Lambeck, [1980 


3.1. Simple tidal dissipation 

We assume that, locally, viscous tidal dissipation scales as the viscosity multiplied 
by the square of the strain rate e. The strain e is the fractional (dimensionless) 
deformation induced by the tide. It is of order of H/R^, where H is the amplitude 
of the tide. The rate of strain is therefore on the order of e = Hf/R^, where / is the 
forcing frequency of the tide. Viscous heating per unit volume scales as pz/e^, which 
for a uniform planet equates to as a whole. The peak-to-peak amplitude of 

the lunar semidiurnal equilibrium tide is 


H 




(15) 


in which M(j is the mass of the Moon, a is the semimajor axis of the Moon’s orbit. 


and h 2 is the amplitude Love number (Kaula, 1968), 


^2 = X 1 + 


19 


/i 


2 pgR® 


-1 


(16) 


The dimensionless quantity {19/2){p/pgR^) compares the material rigidity p (resis¬ 
tance to shear) to the weight of the Earth. The angular frequency of the semidiurnal 
lunar tide — the only tide we will consider here — is / = 2Vt — 2ni, where ni denotes 
the angular velocity of the Moon in its orbit. 
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Viscous heating must be modified to take into account how quickly matter can 
respond to forcing. The resulting heuristic expressionj^ for tidal dissipation in Earth 


IS 




^ (17) 

where Tj, represents the viscous relaxation time scale, which is best written (M. 
Efroimsky, personal communication) in the form 


Tu — Tm + Td 


(18) 


with tm = being the Maxwell time scale and td = l9i'/2gR,^ being Darwin’s 


time scale for relaxation of a viscous globe under its own gravity (Lamb, 1932, p. 641) 


Expressions in the form of Eq. [^are often encountered in the planetary literature to 
describe tidal dissipation in a variety of objects (e.g., 


Moore, 2003 Garrick-Bethell 


et ah, 2010). For the shear modulus /r we assume 


fi = 0 

= 1.3 X 10^^ exp {—a. 
/i = 1.3 X 10^1 


Ti > Tcrit 
T„it >Ti> T„ 
Tso\ > Ti 


(19) 


where the units of /i are kg/m/s. The shear modulus is related to the Love numbers 


^2 (Eq 16) and /c 2 (Eq 23). For Earth today, Kaula (1968) and Lambeck (1980) 


give ^2 ~ 0.6 and k 2 ~ 0.3. These values imply that fi = 1.3 x 10^^ kg/m/s can be 
used as an effective average rigidity of the Earth today. The presumed temperature 


dependence of /i in Eq. 19 is arbitrary. 


Viscous dissipation can be expressed in terms of Q. Kaula (1968, Eq. 2.3.1) 
defines Q as 

97r 1 /• . 

( 20 ) 


2ir 1 / .1, , 


According to Makarov and Efroimsky (2014), this is exact if multiplied by a factor of 19/25. 
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where E represents the energy in the wave and the integral is over one full cycle. 
This is equivalent to writing 

Q = .Epot/i^tide; (21) 

where the total potential energy in the tidal wave is (]Kaula 1968, p. 201). 


Bpo. = \k,GMl ^/, 


where the potential Love number is 


^2 — I 1 + 


19 p 
2 pqRq 


-1 


( 22 ) 


(23) 


We use Eq. to generate 

At temperatures below the solidus, for which viscosity is very large, ordinary 
viscous dissipation at tidal frequencies is probably unimportant. To bring the Moon 
to its present orbit at the present time requires ~ 30, which we obtain by setting 
an upper bound on “viscosity” in Eq. 17 of z/max = 1-3 x 10^^ m^/s. We apply the 


upper bound on v only to tidal dissipation in Eq We do not apply it to the 
Rayleigh number and related quantities pertinent to parameterized convection. 

At high temperatures the ordinary viscosity of mantle liquids is expected to be 
very small. With v = 10““^ m^/s, we estimate that Qe ~ o' 2 T 2 ~ lO^'^, where ct 2 = 
is the relevant natural oscillation frequency of a self-gravitating fluid globe 
(Lamb, 1932, art. 262, Eq. 10) and T 2 = Ri/5z/ is the relevant timescale for viscous 


dissipation in a self-gravitating fluid globe when viscosity is small (Lamb, 1932 


art. 355, Eq. 12). Convecting fluid planets may provide useful guidance. General 
constraints set by tidal evolution of the Galilean satellites suggest that Q of Jupiter is 


between 10® and 10® ( 

Goldreich and Soter, 

1966 

. Astrometric observations suggest 

that Q of Jupiter is ~ 3.5 x 10"^ ( 

Lainey et ah 

2009 

) and Q of Saturn just 1400 


(Lainey et ah, 2012), the latter remarkably low compared to the 7 x 10^ estimated 


by Goldreich and Soter (1966). The uncertainty and variance reminds us that origin 


of dissipation in the giant planets is unknown. Similar arguments have been made 
for binary stars; these suggest that Q is of the order of 10^ — 10® (Meibom and 


[ 
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Mathieu, 2005). On the other hand, at least one hot super Jupiter—the 10 Jupiter 
mass WASP-18b, currently parked in a 0.94 day orbit—may require Q > 10® (Hellier 


et ah, 2009) 


3.2. Simple tidal evolution 

The simplest assumption is that the Moon evolved away from the Earth while 
maintaining a circular orbit. We will make this assumption here. The Moon revolves 
about the Earth at the angular frequency 


n, = yG(A4 + Mj)/o3. 


(24) 


The total Earth-Moon system angular momentum ignoring lunar spin and assuming 
a circular orbit is 

L = J- M(^ y /G {M(^ J- ikf|j ) a, (25) 

in which denotes Earth’s angular rotation rate. The canonical model of the Moon¬ 
forming impact assumes that the angular momentum of the Earth-Moon system has 
been conserved. With L = 0, the Moon recedes at the rate 


a = 




(26) 


With d given by Eq the tidal dissipation is related to the rate of braking of 
Earth’s spin by 


A'tide I 




Q, 


a L-m'-’ P'> 

in which = 0.33 M^R'^ represents Earth’s moment of inertia. Note that is 
negative and d and Etide are positive. 

Figure shows thermal evolution of an Earth after a canonical Moon-forming 
impact that forms a Moon. Other particulars are the same as with the moonless 
impact of Figure Tidal heating extends the duration of the runaway greenhouse 
phase by a few hundred thousand years, and greatly increases geothermal heat flow 
for tens of millions of years afterward, but at hrst glance the differences between 
Figures and are subtle. 
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Figure 3: Thermal evolution of the Earth after a canonical Moon-forming impact. Results are 
shown for the 100 bar BSE-equilibrated atmosphere. The results are very similar to Eigurej^ Tidal 
heating extends the runaway greenhouse phase by a few hundred thousand years compared to the 
moonless case. Geothermal heat flow F is about an order of magnitude higher than in the moonless 
case for a few tens of Myrs. 

Changing the radiative properties of the atmosphere has greater effect than 
adding the Moon. Figure shows a comparable thermal evolution under a 1000 
bar post-impact atmosphere. This is toward the thicker end of the range we might 
expect for Earth, but the thicker atmosphere also serves as a proxy for additional 
opacity sources that are undoubtedly present at very high temperatures but are not 
included in our model, or a proxy for the greater opacity of an atmosphere that is 
more reduced than the BSE. The thicker atmosphere increases how long Earth re¬ 
mains liquid by an order of magnitude, from ~ 0.3 Myrs for the 100 bar atmosphere 
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Figure 4: Thermal evolution of the Earth after a canonical Moon-forming impact, with more 
volatiles. Results are shown for a 1000 bar BSE-equilibrated atmosphere. The thicker atmosphere 
more than the doubles the duration of the runaway greenhouse phase, and increases the Earth’s 
fully molten stage from 0.3 Myrs (Fig. to 3 Myrs. 


to ~ 3 Myrs. The thicker atmosphere also doubles how long the Earth remains in 
the runaway greenhouse state (from ~2 to ~4 Myr). Since the two cases start with 
the same amount of thermal and rotational energy, the differences between Figures 
and 1^ stem from how quickly the atmosphere radiates when it is much hotter than 
the runaway greenhouse state. This in turn is determined by the opacity of the at¬ 
mosphere when the surface is very hot, Tg > 2500 K, a kind of atmosphere that is 
only now beginning to be studied. 

It is interesting to consider what these models predict for the evolution of the 
distance to the Moon and the parameter Qe of the Earth. Figure]^ shows two cases 
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Figure 5: The distance to the Moon a and the tidal dissipation parameter of the Earth after 
a canonical Moon-forming impact. Results are shown for high (dotted) and low (solid) values of 
Q in the fully molten Earth. As it congeals, the mantle becomes extremely dissipative with Qe 
near unity and the Moon’s orbit evolves quickly. After the mantle freezes the viscosity increases, 
Qe increases, and the Moon settles down. Both models assume a relatively thin 100 bar BSE- 
equilibrated atmosphere. The 1000 bar atmosphere gives qualitatively similar results. 


that differ only in the assnmed npper bonnd on 10^ < Qe < 10® for the fnlly liqnid 
Earth. Both cases assnme the Earth is enveloped in the 100 bar atmosphere. There 
are no striking differences between the two models. A minor difference is that the 
higher Q© model delays tidal heating somewhat, thus extending the lifetime of the 
magma ocean (defined as Tj > Tgoi) from 70 to 100 Myrs. The 1000 bar atmosphere 
gives qualitatively similar results, with a somewhat slower tidal evolution and a 
somewhat longer-lived magma ocean. 
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3.3. The Evection and Eviction Resonances 

The evection resonance describes the state in which the perigee of the Moon’s 


orbit processes in one year (Touma and Wisdom, 1998). While the Moon is in the 
evection resonance, its line of apsides remains nearly at a constant angle to the Snn, 
and thns the Snn’s gravitational torqne on the Earth-Moon system does not average 
ont, bnt rather becomes cnmulative. This enables transfer of angnlar momentnm 


from the Earth-Moon system to the Earth’s orbit abont the Snn (Cnk and Stewart 


2012 ). 


The possible importance of the evection resonance to the history to the lunar 


orbit was pointed out by Touma and Wisdom (1998). They showed that, if the 
Earth-Moon angular momentum were conserved, the evection resonance would have 
been encountered at a ~ 4.6/?®, and that capture by the resonance from a nearly 
circular orbit is likely if da/dt < 1 km/yr, but unlikely if da/dt > 10 km/yr. They 
compared these results to a conventional constant Q model that predicts da/dt ^ 150 
km/yr at a ~ 4.6i?0. Touma and Wisdom ( |1998 ) found that if captured the Moon 
can evolve into a highly eccentric orbit with e ~ 0.5, which in their model is what 
allows the Moon to escape the resonance. 

Integrations of the lunar orbit backward in time indicate that today’s 5° inclina¬ 
tion to the ecliptic maps to a ~10° inclination with respect to Earth’s equator when 


the Moon was near the Earth (Goldreich, 1966 Mignard, 1980). In the giant impact 
hypothesis, the Moon should have formed in Earth’s equator and hence should have 


no inclination today. Ward and Canup (2000) proposed that a resonant interaction 
between the nascent Moon and the debris disk interior to it could generate an in¬ 


clination of the right magnitude. Touma and Wisdom (1998) suggested that tidal 
evolution could raise the Moon’s inclination if the Moon were captured successively, 
hrst by the evection resonance, and then by a closely related “eviction” resonance. 
The evection resonance gives the lunar orbit substantial eccentricity. The eviction 
resonance can then give it substantial inclination. The chief challenge is that, to 
be captured by the eviction resonance, the Moon must be evolving inward at the 
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Figure 6: Distance to the Moon a and lunar recession rate da/dt are plotted for our standard case 
with a 100 bar BSE atmosphere. Touma and Wisdom (19981 show that the evection resonance is 
encountered at a ~ 4.6i?0, and that resonance capture from a near circular orbit requires da/dt < 1 
km/yr. The dot-dashed lines indicate that the resonance at 4.6i?0 is met early, after just 20,000 
Myrs (before the mantle has cooled to the liquidus), and that da/dt < 0.3 km/yr, the latter making 
resonance capture very probable. 


time (d < 0), possible only if the Moon were in a very eccentric orbit. They hnd by 
numerical experiment that capture by the eviction resonance requires \da/dt\ < 0.2 
km/yr. 

Figure illustrates that thermal blanketing by Earth’s atmosphere slows early 
tidal evolution of the Earth-Moon system by orders of magnitude compared to 
da/dt = 150 km/yr. The slow rate of tidal evolution makes resonance capture seem 
likely. But we do not see much hope for the eviction mechanism here. The problem 
is that the resonance encounters occur very early, when the Earth is still effectively 
inviscid, but not so early that the airless Moon is likely to be fully molten. With 
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tidal dissipation efficient in the Moon but inefficient in the Earth, the Moon’s orbit 
cannot acquire the substantial eccentricity that it would need to be captured by 
eviction. We expect rather that the Moon does not acquire a substantial eccentricity 
until the Earth begins to freeze and becomes small, which becomes the case when 
20i?® < a < 40i?®. 


3.4- New Conventions 
Recently, 


Cuk and Stewart 


(2012) showed that, if the Earth-Moon system began 
with more angular momentum than it has today, the evection resonance occurs fur¬ 
ther from the Earth (because the faster spinning Earth is more oblate), which makes 
resonance capture more probable. In their model the resonance remains occupied 
during a second contraction phase where both e and a slowly decrease. They also 
showed that, if the Earth-Moon system began with more angular momentum than 
it has today, it would lose that angular momentum while captured in the evection 
resonance. 

To approximate the evection mechanism at a level appropriate for this study, 
we assume that when the Moon encounters the resonance at distance a,., it remains 


trapped while L > Following Touma and Wisdom (1998), we approximate the 
location of the evection resonance by 


Or \ 3 ni 


(28) 


in which ni (Eq. 24) and n 2 are the angular velocities of the month and year, and 

2 


j, - I §- 


(29) 


approximates the oblateness of the Earth with J 20 = 0.00108. With these assump¬ 
tions, 

Etide = - -Md nlordr (30) 

and the angular momentum of the Earth-Moon system decreases as 


L = + Md niOrdr- 


(31) 
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Figure 7: Thermal evolution of the Earth after a high angular momentum Moon-forming impact. 
Results are shown for 100 bar atmospheres. Excess angular momentum in the Earth-Moon system 
is removed by the evection resonance, here crudely modeled by holding the Moon at the resonance 
distance until the angular momentum has been reduced to modern levels. This captures the spirit 
if not the details of Cuk and Stewart’s model. Compare to Figure]^ 


Figure shows thermal evolution of the Earth after a Moon-forming impact. For 
Earth, the liquid magma ocean phase is prolonged to ~ 6 Myr by the added heat 
released by tidal dissipation but otherwise resembles the angular momentum con¬ 
serving case seen in Figure above. 

The evolution of Qe in the Earth is the topic of Figure]^ the analog to Figure]^ 
for a high angular momentum Moon-forming impact. Figure]^ shows lunar recession 
velocities for the middle case from Figure 

In general, tides in the Earth increase the eccentricity of the lunar orbit while 
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Figure 8: Evolution of lunar distance a and Earth’s Q,^ after high angular momentum Moon¬ 
forming impacts. Three cases are shown that correspond to different upper bounds on in the 
fully molten Earth. All cases assume the 100 bar BSE-equilibrated atmosphere. Results can be 
compared to Figure for constant angular momentum. 


tides in the Moon decrease the eccentricity of the lunar orbit (Goldreich and Soter 


1966). Recent discussions of lunar orbital evolution have parameterized this struggle 


as an A factor (Mignard 


1980 


Touma and Wisdom 


1998 


Cuk and Stewart 


2012 ), 


^ ^ 2(1 Q® _ p® R® ^ 2(1 Q(i 


10 


k2(L Q(L 


(32) 


Q(i k2® R(i Q(i k2® Q(i k2® 

that quantihes the relative importance of tidal dissipation in the Moon to tidal 
dissipation in the Earthj^ Today ^ 2(1 = 0.02 and k 2 ® = 0.3 (Kaula, 1968), which 
results in R 0.5 (with = 15). When ^4 !:§> 1, dissipation is relatively strong in 


^The power “5” is incorrectly written “3” 


Touma and Wisdom 

(1998 
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Figure 9: Lunar recession velocities for a high angular momentum Moon-forming impact. The case 
shown is the middle one from Figure]^ In resonance the Moon approaches the Earth {da/dt < 0 
indicated by dotted curve). Earth encounters the evection resonance before it begins to soldify. 
Recession rates are low enough that resonance capture is certain. Details in the structure of da/dt 
reflect changes in the mantle’s rheological properties as it cools. 


the Moon, so that the orbit is circularized. This is expected at early times when Earth 
is fully molten and is very high, while the Moon has little or no atmosphere and 
thus quickly cools to a highly dissipative state with Q(i in the range of 1 to 10. On 
the other hand, when A -C 1, dissipation in Earth is more important than dissipation 
in the Moon, so that tides raise the eccentricity of the lunar orbit. This is the case 
at later times when Earth’s mantle is freezing and Qe 10. 

All models that we have looked at, with or without resonance capture over a 
wide range of atmospheres and planetary albedos, predict (apart from an initial 


transient) the same evolution of A as a function of a (Figure 10 A assumes current 
lunar parameters of /c 2 (i = 0.02 and Q(i = 20). The constancy of A{a) (i.e., of 
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Distance to Moon, a, in Earth Radii 


Figure 10: Evolution of the A factor that parameterizes the relative efficiency of tidal dissipation 
in the Moon to tidal dissipation in the Earth (Eq 32). Lunar parameters fc 2 ([ = 0.02 and Q(^ = 20 
are held fixed. Several cases are shown: the conventional case from Fig. (dashed line), the three 
“evection” cases from[^ and another with a 1000 bar atmosphere (solid lines). Differences at the 
beginning stem from different assumptions about in the fully molten Earth. Otherwise the 
course of A(a) is the same in all our models. In all cases the best prospects for growing a lunar 
eccentricity is with A^\, when the distance to the Moon was between 20 and 40 i?®. 


as ^ 2 © does not stray from 1.5 until late) implies that the tidal evolution of Q©(a) is 
determined by a general principle, probably conservation of energy, rather than by 
details of our model. 

A second key point made in Figure [T^ is that, in all cases, A ;:§> 1 when the 
Moon encounters the evection resonance. In the class of models being considered 
here, there appears to be no way for tidal evolution to be both slow enough to hold 
the Moon within a < 7R® yet simultaneously to be highly dissipative (low Q©). The 
reason for this is that small values of Q© cause fast tidal evolution, quickly raising a 
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above 15i?0. Thus the Moon could not acquire much eccenticity while in the evection 
resonance. 


A third point made in Figure 10 is that, in all cases, A -C 1 when 20<a< 
4Oi?0. In this phase of lunar orbital evolution dissipation in the Earth would have 
been much more important than dissipation in the Moon, so that tides raised in the 
Earth by the Moon would have raised the eccentricity of the lunar orbit. How big 


e can grow is unclear. Touma and Wisdom (1998) found that e = 0.5 is achievable 
in the evection resonance. Presumably tidal flexure in the Moon caused by a highly 
eccentric orbit would at some point trigger considerable tidal melting within the 
Moon and a consequent drop of Q(i. Later still, after the Earth has mostly solidihed, 
A must rise again and most of any remaining eccentricity is taken out of the Moon’s 
orbit. 

How a late episode of lunar eccentricity might be twisted into today’s lunar 


inclination is another puzzle beyond the scope of this study. Touma (2000) mentions 
a secular resonance for growing inclination when the Moon was at 28i?0 if e > 0.6, 
but the mechanism is not otherwise described and we can say no more. We do 
note that, in general, prospects for giving the Moon an inclination are probably best 
when the Moon’s orbit was neither wholly under the Sun’s influence nor the Earth’s 
influence, which roughly corresponds to 15i?0 < a < 3Oi?0. During this interval 
the Laplace plane coincides with neither the ecliptic nor Earth’s equator which, if 


nothing else, creates an opportunity for helpful accidents (Ward, 1981). 


4. Conclusions 

The greenhouse effect of a thick water-rich atmosphere slows the rate that a 
planet cools after a giant impact. Here we show that for water-rich atmospheres, 
the difference between the runaway greenhouse limit of ~ 280 W/m^ and absorbed 
sunlight (which for the faint young Sun and and a planetary albedo of 0.3 is ~ 170 
W/m^) sets an upper bound of ~ 110 W/m^ on how quickly the Earth can cool. 
This in turn sets a 10® year timescale for a magma surface after an impact on the 












scale of the event that formed the Moon. 

Tidal heating is a major term in the early Earth’s energy budget. Tidal heating 
occurs mostly in mantle materials that are just beginning to freeze, which frustrates 
freezing. The atmosphere’s control over the planet’s cooling rate sets up a negative 
feedback between viscosity-dependent tidal heating and the temperature-dependent 
viscosity of the magma ocean. While this feedback holds, the rate that the Moon’s 
orbit evolves is limited by the modest radiative cooling rate of Earth’s atmosphere, 
which in effect tethers the Moon to the Earth. Consequently the Moon’s orbit 
evolves orders of magnitude more slowly than in conventional models that describe 


tidal evolution through hxed values of Q (Touma and Wisdom 


1998 


Cuk and Stew¬ 


art, 2012). The slow tidal evolution ensures that the evolving Earth-Moon system 


will be captured by the strong evection resonance, but in our models the resonance 
encounter always occurs when tidal dissipation is more efficient in the Moon than 
in the Earth, which suggests that capture by the evection resonance is unlikely to 
generate much eccentricity in the lunar orbit. Whether signihcant angular momen¬ 
tum can be extracted from the Earth-Moon system by cumulative solar torques from 
a low eccentricity lunar orbit is a question beyond the scope of this research, but 
the answer to this question would have implications for Cuk and Stewart (2012)’s 
inventive hypothesis. 

Slow orbital evolution also promotes capture by other, weaker, higher order res¬ 
onances (Touma and Wisdom, [1998 ; Touma, 2000). Resonance capture by one or 
more of these resonances may have play a key role in the orbital evolution of the 


Earth-Moon system, a possibility raised but not elaborated upon by Touma and 


Wisdom (1998). Tidal dissipation in the Earth is most efficient at later times when 


the Earth’s mantle has cooled enough to develop signihcant viscosity. In all our 
models, tidal dissipation is most efficient in the Earth when the Moon is between 
2Oi?0 and 4Oi?0. During this interval it is likely that the Moon’s eccentricity grew 
considerably. The upper limit on eccentricity may be set by vigorous tidal heating in 
the Moon — once the Moon melts, tidal dissipation in the Earth and the Moon will 
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be more comparable. We speculate that this period of high eccentricity set the table 
for some other process, perhaps capture by an otherwise unremarkable, little-noted 
resonance, that exploited its opportunity by leaving its mark on the inclination of 
the Moon. 
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